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ABSTRACT 

There are two kinds of expressions on the precession of orbital plane of a binary 
pulsar system, which are given by Barker & O'Connell (1975) and Apostolatos et al. 
(1994), Kidder (1995) respectively. This paper points out that these two kinds of or- 
bital precession velocities are actually obtained by the same Lagrangian under different 
degrees of freedom. Correspondingly the former expression is not consistent with the 
conservation of the total angular momentum vector; whereas the latter one is. Damour 
Sz Schafer (1988) and Wex & Kopeikin (1999) have applied Barker & O'Connell's orbital 
precession velocity in pulsar timing measurement. This paper applies Apostolatos et 
al. & Kidder's orbital precession velocity in pulsar timing measurement. We analyze 
that Damour & Schafer's treatment corresponds to negligible Spin-Orbit induced pre- 
cession of periastron. Whereas the effects corresponding to Wex & Kopeikin and this 
paper are both significant (however they are not equivalent). The observational data of 
two typical binary pulsars, PSR J2051-0827 and PSR J1713-I-0747 apparently support 
significant Spin-Orbit coupling effect. Further more, the discrepancies between Wex 
&: Kopeikin and this paper can be tested on specific binary pulsars with orbital plane 
nearly edge on. If the orbital period derivative of double-pulsar system PSRs J0737- 
3039 A and B, with orbital inclination angle i = 87.7^29deg, is much larger than that of 
the gravitational radiation induced one, then the expression of this paper is supported, 
otherwise Wex &: Kopeikin's expression is supported. 

Subject headings: pulsars: binary pulsars, geodetic precession: individual (PSR J2051- 
0827, PSR J1713+0747, PSRs J0737-3039 A and B) 



1. Introduction 

In the gravitational two-body problem with spin, each body is precessing in the gravitational 
field of its companion, with precession velocity of 1 Post-Newtonian order (PN) (Barker &: O'Connell 
1975, hereafter BO). This precession velocity is widely accepted. But how the orbital plane reacts 
to the torque caused by the precession of the two bodies has two kinds of treatments. BO's orbital 
precession velocity is obtained by assuming that the angular momentum vector, L, precesses at the 
same velocity as the Runge-Lenz vector, A. 
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On the other hand, in the study of the modulation of the gravitational wave by Spin-Orbit 

(S-L) coupling effect on merging binaries (Apostolatos ct al. 1994, Kidder 1995, hereafter AK) 
obtain an orbital precession velocity that satisfies the conservation of the total angular momentum, 
J, and the triangle constraint, J = L + S (where S is the sum of spin angular momenta of two 
bodies. Si and S2). 

The discrepancy in the assumptions leads to different behaviors in physics between BO and 
AK's orbital precession velocities. The former isn't consistent with the triangle constraint, whereas 
the latter is. And the reason of this discrepancy is that the former actually assumes that the 
four vectors. Si, S2, ri and r2 (ri, r2 are position vectors of the two bodies respectively), are 
independent; whereas, the latter assumes that the independent vectors are either Si, ri and r2; or 
S2, ri and r2. 

The discrepancy between BO and AK's precession velocities is analogous to the following 
case. The motion of a small mass at the bottom of a clock pendulum can be described m x — y 
plane. However if we treat the degree of freedom of this small mass as 2, then this small mass 
can move every in the 2-dimension space, and the length of the pendulum is not a constant. In 
other words, once the length of the pendulum is a constant, the degree of freedom is 1 instead of 
2. Correspondingly if the free vectors of a binary system is 4, as treated by BO, then the triangle 
constraint cannot be satisfied (or J cannot be a constant vector). On the other hand, if the triangle 
constraint is satisfied, the number of free vectors is 3, instead of 4. 

In the application to pulsar timing measurement, BO's orbital precession velocity is treated 
in two different ways, by Damour & Schafer (1988) and Wcx & Kopcikin (1999) respectively. 
The former one predicts insignificant S-L coupling induced precession of periastron, co^ , and the 
derivative of orbital period, P^; whereas the latter predicts significant u>^ and P^. The paper points 
out that the discrepancy is due to Damour &; Schafer and Wex &; Kopeikin calculated effects in 
different kinds of coordinate systems, the former coordinate system is not at rest to "an observer" 
at the Solar System Baryon center (SSB), whereas, the latter is at rest to it. In other words, 
the former coordinate system has non-zero acceleration to SSB; whereas, the latter one has zero 
acceleration to SSB. 

This paper calculates observational effect corresponding to AK's orbital precession velocity, 
which uses the same coordinate system as that of Wex & Kopeikin (1999). Significant u>'^ and Ph 
are given, but they are not equivalent to the results given by Wex & Kopeikin (1999). The validity 
of these two expressions can be tested by specific binary pulsars with orbital inclination close to 
7r/2, i.e., the double-pulsar system PSRs J0737-3039 A and B. 

This paper contains four parts: (a) the physical discrepancy between BO and AK's orbital 
precession velocity (Sect 2,3); (b) the derivation of S-L coupling induced effect corresponding to 
AK's expression of orbital precession (Sect 4,5,9); (c) the discrepancy between the coordinate 
systems used by Damour & Schafer (1988) and Wex &: Kopeikin (1999), as well as the relationships 
among three kinds of S-L coupling induced effects, Damour & Schafer and Wex &; Kopeikin and 
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this paper (Sect 6); (d) confrontation of the three S-L couphng effects with observational data of 
PSR J2051-0827 and PSR J1713+0747, and different predictions on uj^ and of PSRs J0737- 
3039 A and B by different S-L couphng models (Sect 7,8). 



2. Orbital precession velocity 

This section introduces the derivation of the orbital precession velocity of BO and AK. 



2.1. Derivation of BO's orbital precession velocity 

BO's two-body equation was the first gravitational two-body equation with spin, which consists 
of two parts, the precession velocity of the spin angular momenta vectors of body one and body 
two, and the precession velocity of the orbital angular momentum vector. Body one precesses in 
the gravitational field of body two, with precession velocity (BO), 



. L(4 + 3m2/mi) , ^2 
"i = 



S2-3(L-S2)L (1) 



where L, Si and S2 are unit vectors of the orbital angular momentum, the spin angular momentum 
of star 1 and star 2, respectively. 

tl2 can be obtained by exchanging the subscript 1 and 2 at the right side of Eq.(l). The first 
term of Eq.(l) represents the geodetic (de Sitter) precession, which corresponds to the precession 
of Si around L, it is IPN due to ^ ~ (^)^(^); and the second term represents the Lcnse-Thirring 
precession. Si around S2, which is ^ times smaller than the first, therefore, it corresponds to 1.5PN. 
The the precession velocity of the spin angular momenta vectors is confirmed by other authors using 
different method. 

However for the precession velocity of the orbit, there are different expressions. BO's orbital 
precession velocity is given as follows. 

The total Hamiltonian for the gravitational two-body problem with spin is given (BO, Damour 
& Schafer 1988), 

H = Hn + HiPN + H2PN + Hs , (2) 

where H^, Hipn and H2PN are the Newtonian, the first and second order post- Newtonian terms 
respectively. Hs is the spin-orbit interaction Hamiltonian (BO, Damour &; Schafer 1988), 

il5 = E(2 + 3^)(^), (3) 

a=l 

where a -|- 1 is meant modulo 2 (2-|-l=l), mi, m2 are the masses of the two stars, respectively, 
r = a(l — e^)^/^, a is the semi-major axis, e is the eccentricity of the orbit. Notice we uses G = c = 1 
until discussing observational effects in Sec 5 to Sec 7. 
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The BO equation describes the secular effect of the orbital plane by a rotational velocity vector, 

0,s, acting on some instantaneous Newtonian ellipse. Damour & Schafer (1988) computed Qs in 
a simple manner by making full use of the Hamiltonian method. The functions of the canonically 
conjugate phase space variables r and p are defined as 

L(r, p) = r X p , (4) 

A(r, p) = p X L - GMn"^- , (5) 

r 

where r = ri — r2, M = mi + m2, n = mim2/M. The vector A is the Runge-Lenz vector (first 
discovered by Lagrange) . The instantaneous Newtonian ellipse evolves according to the fundamental 
equations of Hamiltonian dynamics (Damour k, Schafer 1988) 

t = {L,H], (6) 

A = {A,i7}, (7) 

where {, } denote the Poisson bracket. L and A are first integrals of i^^v, only Hipn + H2PN + Hs 
contributes to the right-hand sides of Eq.(6) and Eq.(7), in which Hip^ + H2PN determines the 
precession of periastron, in IPN it is given, 

" - na(l - e2) ' 

where Ph is the orbital period. To study the spin-orbit interaction, it is sufficient to consider Hs- 
Thus replacing H of Eq.(6) and Eq.(7) by Hs obtains (Damour & Schafer 1988), 

{^)s = {L,Hs} = n*sSxL, (9) 
flA . . . ^ . 

{—)s = {A, Hs} = n*s[S - 3(L . S)L] x A . (10) 

where 

5(4 + 3m2/mi) 
^~ 2^ ■ 

By Damour &: Schafer (1988), S represents a linear combination of Si and S2. For simplicity of 
discussion, and for consistence with Wex & Kopeikin's application of fls (Wex & Kopeikin 1999), 
we assume S = Si (the other spin angular momentum is ignored) until Sec 4 where the general 
binary pulsar is discussed. 

The solution of Eq.(9) and Eq.(lO) gives the S-L coupling induced orbital precession velocity 
(Damour & Schafer 1988) 

Qs = ^*s[^ - 3(L • S)L] . (12) 
By Eq.(9) and Eq.(12), the first derivative of L can be obtained, 

^ = n*sSxt, (13) 
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and by Eq.(l) the first derivative of S (recall S = Si) can be written 

^ = ji^L X S = n*s^L X S , (14) 

where Cl^ is the first term at the right hand side of Eq.(l). By Eq.(13) and Eq.(14) L precesses 
slowly around S, 1.5PN, as shown by Eq.(13); whereas S precesses rapidly around L, IPN, as 
shown in Eq.(14). Therefore the BO equation predicts such a scenario that the two vectors, L 
and S precess around each with very different precession velocities (typically one is larger than the 
other by 3 to 4 orders of magnitude for a general binary pulsar system). 



2.2. Orbital precession velocity in the calculation of gravitational waves 

In the study of the modulation of precession of orbital plane to gravitational waves, the orbital 
precession velocity is obtained in a different manner and the result is very different from that given 
by Eq.(12). 

Since the gravitational wave corresponding to 2.5PN, which is negligible compared with S-L 
coupling effect that corresponds to IPN and 1.5PN, the total angular momentum can be treated 
as conserved, J = 0. Then the following equation can be obtained (BO), 

J7o X L = -fJi X Si - J72 X S2 . (15) 

Notice that as defined by BO and AK, L = /xM^/^r^/^L. In the one-spin case the right hand side 
of Eq.(15) can be given (Kidder 1995), 

and considering that J = L + S, Eq.(16) can be written, 

S = -L(4 + ^)(JxS) . (17) 
2r-^ mi 

Prom Eq.(15), L can be given. 

By Eq.(17) and Eq.(18), L and S precess about the fixed vector J at the same rate with a precession 
frequency approximately (AK) 

"0 = ^3(4 + ^) ■ (19) 
2r^ mi 

Eq.(19) indicates that in the 1 PN approximation, L and S can precess around J rapidly (IPN) 
with exactly the same velocity. Notice that the misalignment angles between L and S {Xls), L 
and J (Alj) are very different, due to S/L <^ 1 , Xlj is much smaller than A^^. 

Thus, AK's equations, Eq.(18) and Eq.(17), correspond to a very different scenario of motion 
of S, L and J from that given by BO equation shown in Eq.(13) and Eq.(14). 
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3. Physical differences between BO and AK 

This section compares two different scenarios corresponding to BO and AK's orbital precession 
velocity, and points out that BO's orbital precession velocity is actually inconsistent with the 
definition of the total angular momentum of a binary system. 

Section 2 indicates that BO and AK derived the orbital precession velocity in different ways, 
therefore two different orbital precession velocity vectors are obtained, as shown in Eq.(12) and 
Eq.(19), respectively, which in turn correspond to different scenarios of motion of the three vectors. 
This section analyzes that the discrepancy between BO and AK is not just a discrepancy corre- 
sponding to different coordinate systems. Actually there is significant physical differences between 
BO and AK. The total angular momentum of a binary system is defined as BO 

J = L + S , (20) 

Eq.(20) means that J, L and S form a triangle, and therefore, it guarantees that the three vectors 
must be in one plane at any moment. For a general radio binary pulsar system, the total angular 
momentum of this system is conserved in IPN. Therefore we have 

j = . (21) 

Eq.(21) means that J is a constant during the motion of a binary system. Eq.(21) and Eq.(20) 
together provide a scenario that the triangle formed by L, S and J determines a plane, and the 
plane rotates around a fixed axis, J, with velocity Oq- This scenario is shown in Fig 1. which can 
also be represented as 

j = Oo J X L + Oo J X S = , (22) 

Smarr k, Blandford (1976) mentioned the scenario that L and S must be at opposite side of J 
at any instant. Hamilton and Sarazin (1982) also study the scenario and state that L precesses 
rapidly around J. Obviously the orbital precession velocity given by Eq.(19) can satisfy the two 
constraints, Eq.(20) and Eq.(21) simultaneously. 

Can the BO's orbital precession velocity given by Eq.(12) satisfy the two constraints, Eq.(20) 
and Eq.(21) simultaneously ? From Eq.(12), Eq.(13) and Eq.(14), the first derivative of J can be 
written (BO) 

j = ^5 X L + fJ^ X S = J7|S X L + ii^L X S = , (23) 

and since Eq.(20) is defined in BO's equation, then it seems that the BO equation can satisfy both 
Eq.(20) and Eq.(21). 

But in BO's derivation of ttg (Eq.(6)~Eq.(12)), Eq.(20) is never used. The corresponding tls 
can make J = L + S = 0, as shown in Eq.(23), however it cannot guarantee that J = L + S is 
satisfied. In other words, when J ^ L + S, Eq.(23) is still correct. This can be easily tested by 
putting L' = L + aS, or S' = S + /3L (a and (3 are arbitrary constants) into Eq.(23) to replace L 
and S, respectively, obviously in such cases, Eq.(23) is still satisfied (J = 0). 
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Contrarily in AK's derivation of i7o(Eq.(16)-Eq.(18)), the relation Eq.(20) is used. And if we 
do the same replacement of L' = L + aS, or S' = S + /3L in Eq.(22), then Eq.(22) is violated (J 7^ 0). 
This means that for AK's fto, if Eq.(20) is violated then Eq.(21) is violated also. Thus in AK's 
expression, the conservation of the total angular momentum is dependent on Eq.(20), whereas, in 
BO's expression, the conservation of the total angular momentum is independent of Eq.(20). If we 
rewrite Eq.(20) as, 

J = L + S + C , (24) 

then in BO's expression, the conservation of the total angular momentum can be satisfied in the 
case that C 7^ in Eq.(24); whereas, in AK's expression, the conservation of the total angular 
momentum is satisfied only when C = in Eq.(24). This means that the discrepancy between 
BO and AK's orbital precession velocity is physical. It is not just different expression in different 
coordinate systems or relative to different directions. 

Moreover Eq.(22) and Eq.(19) correspond to the following orbital precession velocity, 

no = n*s{s + . (25) 

Obviously Eq.(25) is not consistent with BO's Eq.(12), which demands that the coefficient of the 
component along L be 7 = — 3(L • S), instead of 7 = ^ as given by Eq.(25). 

In other words, once Eq.(20) is satisfied, BO's orbital precession velocity of Eq.(12) must be 
violated. Therefore, BO's orbital precession velocity cannot be consistent with BO's definition, 
J = L + S. 

Actually Eq.(25) can be consistent with Eq.(9), however it is contradictory to Eq.(lO). The 
reason of introducing Eq.(lO) is that without it, Eq.(9) alone cannot determine a unique solution. 

Whereas, Eq.(22) and Eq.(19) can be regarded as solving this problem by using Eq.(9) and 
Eq.(20) instead of Eq.(9) and Eq.(lO) to obtain the orbital precession velocity. 

As defined in Eq.(4) and Eq.(5), L and A are vectors that are determined by different elements 
in celestial mechanics, L(ri,i) and A(Q,i,uj,e) respectively. And these two vectors satisfy different 
physical constraints, i.e., L satisfies Eq.(20) and Eq.(21), whereas, A doesn't satisfy these two 
constraints. 

Therefore, it is conceivable that L and A should correspond to different precession velocities, 
as given by Eq.(9) and Eq.(lO), respectively. However since the discrepancy is only in the L 
component, which does not influence the satisfaction of the conservation equation, Eq.(23), thus the 
discrepancy seems unimportant. And therefore, the precession velocity of L is treated equivalently 
to that of A's, thus the components in L are both treated as 7 = — 3(L • S). Whereas, as given by 
Eq.(25), L component must be 7 = ^ if the triangle constraint is to be satisfied. Therefore, the 
violation of the triangle constraint is inevitable under the assumption that L and A precess at the 
same velocity. 
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4. S-L coupling induced effects derived under different degree of freedom 

As analyzed in Sect 2 and Sect 3, whether the triangle constraint is satisfied or not results 
discrepancy in the orbital precession velocity. This section further analyzes that it is the discrepancy 
in degree of freedom used by BO and AK that leads to the violation or satisfaction of the triangle 
constraint. 

To discuss the S-L coupling induced effects on observational parameters, one need to obtain 
the variation of the six orbital elements under S-L coupling. The way of doing this is from the 
Hamiltonian (corresponding to S-L coupling) to equation of motion, and then through perturbation 
methods in celestial mechanics to obtain S-L coupling induced variation of the six orbital elements, 
and finally transform effects to observer's coordinate system. In this section this process is per- 
formed in the case that the free vectors of a binary system (with two spins) is 3, in which the 
triangle constraint is satisfied. 

It is convenient to study the motion of a binary system in such a coordinate system (J- 
coordinate system), in which the total angular momentum, J is along the z-axes and the invariance 
plane is in the x-y plane. The J-coordinate system has two advantages. 

(a) Once a binary pulsar system is given, A^j, the misalignment angle between J and L, can 
be estimated, from which Jl and u can be obtained easily in the J-coordinate system, which are 
intrinsic to a binary pulsar system. 

(b) Moreover, the J-coordinate system is static relative to the line of sight (after counting out 
the proper motion). Therefore transforming parameters obtained in the J-coordinates system to 
observer's coordinate system, S-L coupling induced effects can be obtained reliably. 



Prom Eq.(3), the S-L coupling induced Hs contains potential part only, therefore we have 
Hs = U, where 

1 /^o , 3m2 g _^ 3mi , 
2mi 2m2 

which can be written as 



= C/l + C/2 = -3 (2S + ;^Sl + ;^S2) • L , (26) 



C/ = ^((TiS-L + a2S2-L) , (27) 
where 

ai=2 + - — , 0-2 = 2 + -( ) . (28) 

2 mi 2 m2 mi 

From Eq.(27) we can have the Lagrange corresponds to S-L coupling, 9 = — [/. And the Lagrange 
equation is 

|(|^)-|^ = ' ('^ = 1'2, ...,/?) (29) 

where is the generalized coordinate (/? is the number of degrees of freedom), given by r^"\ i^2^\ 
Si^\ 52°^ (which represent the position of body 1 and body 2; the spin angular momentum of body 
1 and body 2 respectively, a = 1, 2, 3). 
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Since the times scale of spin down of a pulsar in a binary pulsar system is much longer than 
that of the period of geodetic precession, the magnitude of Si and S2 can be treated as constant 
(these two vectors can only vary in direction). However since the misalignment angle between 
Si and L (Alsi), as well as S2 and L {Xls2), arc also constants in BO's gravitational two-body 
equation (Barker & O'Connell 1975), we have ^(Si • L)/dS[''^ = and 8(82 ■ L)/95^"^ = 0. Thus 
we have d^/dq^, = for = -S^^^ £'2"^ And since Si and S2 are not appeared in the Lagrange S5, 
we have d{d'^ / dq^) / dt = for q^ = s["\s^''\ 

Therefore, we only need to calculate d{d'ss / dq^) / dt and d^/dq^, of Eq.(29) in the case that 
q^ = rj^^^rg"^ . The first term at the left hand side of Eq.(29) corresponds to generalized force, 
which can be written d{d'^ / dq^) / dt = F = j-iaso] and the second term is d^/dq^ = —VU (where 
V represents gradient). Thus Eq.(29) can be rewritten, 

a.. = --Vt/ = -i[aiV^ + a2V^^] , (30) 

The triangle constraint given by Eq.(20) indicates that S and L are dependent. Therefore, we 
cannot treat the free variables of Eq. (30) as rj""* , r^^'' , and . 

Classical mechanics tells us that constraints actually reduce the number of degrees of freedom 
of a dynamic system. The geometric constraint J = L + S can be imposed in Eq. (30) through the 
replacement, S = J — L. Thus the free variables in Eq.(30) is either r^^\ rg"^ and -Sl*^^; or 
r^2^\ and <S'2"^ (depending on different definition of ai and 02 in Eq.(28)). Therefore, performing 
the replacement S = J — L in Eq.(30) actually means that the triangle constraint is imposed on 
the motion of binary system, and the degrees of freedom are reduced from 12 (4 x 3) to 9 (3 x 3). 

Considering that Si and S2 can vary in direction only, a of 5*1°^ S^"^ is given a = 1, 2. Thus 
the degree of freedom of a binary pulsar system without the triangle constraint is 10, and with the 
constraint is 8. 

This is analogy to the calculation of equation of motion of a simple clock pendulum. The motion 
of a small mass at the bottom of a clock pendulum can be described in x — y plane. However if 
we treat the degrees of freedom of this small mass as 2, then this small mass can move every in 
the 2-dimension space, and the length of the pendulum is not a constant. In other words, once the 
length of the pendulum is a constant, the degree of freedom is 1 instead of 2. Correspondingly if the 
degree of freedom of a binary system is 10, then the satisfaction of the triangle constraint cannot 
be guaranteed (or J=const vector cannot be guaranteed). Contrarily, if the triangle constraint is 
satisfied, the degree of freedom 8 instead of 10. By the replacement, S = J — L, Eq.(30) can be 
re- written 

^(S-L) _^[(J-L).L] ^(J.L) ^(L-L) 



1^3 y^3 y^3 ^3 



(31) 



where 



^y^_3(L^_3(Jxr)(^^^^(J^ (32) 

y^O .^0 lytJ 

„(L-L) (L-L)r „(LxV) 

rj^tj rjT'tJ ij^tj 



-10- 



where V is the velocity of the reduced mass. Put Eq.(32) and Eq.(33) into Eq.(31) and finally into 
Eq.(30), we have, 

3 

a^o = ^ki(J - L) • (n X V)n + (72S2 • (n x V)n] 

+^[2cJi{Y x J) - c7i(V X (J - L)) + a2(V X S2)] 

+^^73^[^i(J ^ + ^2(82 X n)] , (34) 
where n is the unit vector of r. Replacing J — L by S, Eq.(34) can be written, 

3 

Blso = ^[o"iS • (n X V)n + 0-282 • (n x V)n] 

+ ^[2(71 (V X J) - (71 (V X 8) + C72(V X 82)] 

+^^^3^[ai(J X n) + (72(82 X n)] . (35) 

If one calculates a^o directly by Eq.(30) without imposing the triangle constraint, then the corre- 
sponding result can be given by replacing J of Eq.(35) by 8, 

3 

^so = ^kiS • (n X V)n + (7282 • (n x V)n] 



+^[(7i(Vx8) + a2(Vx82)] 

+ ^''^3""^ [(71 (8 X n) + (72(82 X n)] . (36) 

Notice that Eq.(36) is equivalent to the sum of Eq.(52) and Eq.(53) given by the BO equation. The 
difference between Eq.(35) and Eq.(36) indicates that whether the triangle constraint is satisfied 
or not can lead to significant differences in a^o, which in turn results in significant differences on 
the predictions of observational effects as discussed in the next section. 

Having a^o, we can use the standard method in celestial mechanics to calculate, 

^ = a,o-n, T = a,„-t, 1^ = a,o • L, (37) 

from which one can calculate the derivative of the six orbit elements and then transform to the 
observer's coordinate system to compare with observation. The unit vector, n, in Eq.(35) to Eq.(37) 
is given by, 

n = Pcos/ + Qsin/ , (38) 
and t is the unit vector that is perpendicular to n, 

t = -Psin/ + Qcos/ , (39) 
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and V = p//x, is given by 

V = --Pcos/ + -Q(e + cos/) , (40) 
P P 

where / is the true anomaly, p is the semilatus rectum, p = a{l — e^), and h is the integral of area, 
h = r^f. P is given by three components, 

= cosr^cosu; — sinrisinwcos Alj , 

= sinJlcoso; + cosfisinwcos Alj , 

P2 = sina;sin Aij , (41) 

and Q is given by three components, 

Qx = — cos sin a; — sin cos u cos A lj , 

= — sinJlsino; + cosOcosa;cos Alj , 

= cos a; sin Al J , (42) 
The unit vector of L and (k = 1, 2) are given, 

L = (sinALjcosr/L, sinAijsinryi, cosAlj)^ , (43) 

S«; = (sinAj5«cosr75«, sinAj5«;Sin?75K, cosAj5k)^ . (44) 

In the perturbation equation, the acceleration of Eq.(35), ago, is expressed along n, t and L respec- 
tively. We can use ai, 3.2 and as to represent terms corresponding to the three terms containing 
brackets [, ] at the right hand side of Eq.(35), respectively. Projecting ai onto L, we have 

Wi = a.i-t = ^[SxiriyVz - n^Vy) + Sy{nzVx - rixVz) + Sz{nzVy - %V;)] 

{hx sin Ai J cos riL + Uy sin Xlj sin riL + Uz cos Xlj) , (45) 

where Ux, Uy, Uz and Vx, Vy, Vz are components of n and V along axes, x, y and z, respectively. 
Similarly, projecting a2 onto L, we have, 

= a2 • L = ^(V'j/JsinALjcosr/L - T^-^sinALjsinryi,) 

+[^ {VxSy - VySx) + {VxS2y - VyS2x)] COS Xlj . (46) 
Finally as can also be projected onto L, 

= as • L = ^^^[a^{v;Sx - V^Sy) + a2{V;S2x - 1^;^2,)] , (47) 
where = 3rna;, Vy = 3rny and r = ^ sin/. Therefore, the sum of W is 

W = Wi + W2 + W3 . (48) 
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The effect around J can be obtained by perturbation equations (Roy 1991, Yi 1993, Liu 1993) and 
Eq.(48) _ 

dn _ Wr sin(u; + /) 1 
dt na^Vl - sinA^j 
where n is the angular velocity. Averaging over one orbital period we have 

< ~^ >= r. -^/-i ^^Ta^/?"^- X {PzSmU + cos Uj)[{PyQz - PzQy){S:a(yi + S2x(T2) 

dt za'^^l — e^)'^/^ sm Alj 

+{P,Qx - PxQz){Syai + S2ya2) + {P,Qy - PyQ:c){S,ai + 82,02)] ■ (50) 

Notice that the average value of Eq.(50) depends on Wi only, the contribution of W2 and VF3 to it 
is zero. With S/ sinA^j ^ L, we have dil/dt ^ L/a^, which corresponds to IPN. 

The duj/dt Ccin b6 obtain6ci by Ccilciilcition of S — ^so * 

n and T = a^o • t. Since ai • n and ai • t 
are 1.5PN. Therefore, it is sufficient to consider the projection of a.2 and as onto n, t respectively, 
thus we have, 

<(a2-n)cos/>=^^^^^i^^(P,Q,-P,g,) , (51) 

< (a, • t) sin / >= :%/2/ j(P^Qy - PyQ^) ' (52) 

<(a3-t)sin/>=^^^^^l^^(P.Q,-P,Q.) , (53) 

< (as • n) cos / >= . (54) 

From Eq.(51) to Eq.(54), we have 



d^ = ^/LZ{[_a2 . n] cos / + (1 + -)[a2 • t] sin/} 
dt nae p 

= 2(1 - e2)3/2a3 ^Py^"" ~ P^^y^ ■ (55) 

therefore, by the standard perturbation(Roy 1991, Yi 1993, Liu 1993), the advance of precession of 
periastron induced by S-L coupling is given by 

duj duj' dVL 

By putting Eq.(51) and Eq.(52) into Eq.(56), and averaging over one orbital period we have, 

duj laiJ r> ^ X 

< -dt 2(l-e^)W ^'^^^- - - -dt • 

Using perturbation equations as in (Roy 1991, Yi 1993, Liu 1993), and by Eq.(49) and Eq.(57), we 
have 

d'cu 7aiJ ^ „ ^ N d^} 9 Arj 

< Ht >= 2(i-e^)3/v <^'g' - + 'it ^ 
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- 2(1 - e2)3/2«3-(^^Q- - P^Qy) + Oic-') . (58) 

Eq.(50), Eq.(57) and Eq.(58) indicate that the magnitude of dQ,/dt, and dw/dt are both L/a^ 
(IPN), whereas, duj/dt can be 1.5PN (or zero in IPN) by Eq.(57). 

dn/dt (IPN) of Eq.(49) is equivalent to $s (IPN) which is given by Wex & Kopeikin (1999). 
This is because the averaged value of dQ/dt depends only on ai, the first term containing bracket 
[,] in ago, as shown in Eq.(35). And both this paper (a^o of Eq.(35)) and that of BO equation (a'^^ 
of Eq.(36)) give the same ai. Thus, different authors give the equivalent value on the averaged 
dn/dt. 

Whereas, duj/dt of Eq.(57) and 4>s (IPN) given by Wex & Kopeikin (1999) are very different in 
magnitude. The difference is due to the fact that dw/dt given by Eq.(57) of this paper is obtained 
by the a^o of Eq.(35); whereas, the corresponding duj/dt of Wex and Kopeikin (1999) is obtained 
by the a^^ of Eq. (36) , which is equivalent of replacing J of Eq. (35) by S. 

And in turn, the difference between ago and a'gg is due to that ago satisfies the triangle con- 
straint; whereas a'^^ doesn't. Therefore, a small difference in the equation of motion causes signifi- 
cant discrepancy in the variation of elements, such as dw/dt. 



5. Effects on Co, x, Pf, 

As shown in Sect 4 the S-L coupling effect can be treated as a perturbation to the Newtonian 
two-body problem, and by the standard method in celestial mechanics the variation of six orbital 
elements can be obtained. This section calculates what w, x, Pt, are for an observer when the 
variation of the six orbital elements is given. The results of this section is independent of the 
degree of freedom used in a binary system. 

The observational effect is studies in such a coordinate system that the vector Kq (corre- 
sponding to line of sight) is along the z'-axis; the x' axis is along the intersection of the plane of 

the sky and the invariance plane; and the y'-axis is perpendicular to x' — z' plane, as shown in 
Fig 2a. Obviously this coordinate system is at rest to "an observer" at the SSB. The relationship 
of dynamical longitude of the ascending node, il, dynamical longitude of the periastron, co, and the 
orbital inclination, i, can be given (Smarr & Blandford 1976, Wex Sz Kopeikin 1999) 

cos i = cos / cos Xlj — sin A lj sin / cos Q , (59) 

and 

sin i sin uj°^^ = (cos / sin Xlj + cos A lj sin / cos O) sin a; 
+ sin / sin J7 cos a; , 

sin i cos lo°^^ = (cos / sin Xlj + cos A lj sin / cos cos w 



(60) 
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— sin 7 sin sin a; , (61) 

where I is the misahgnment angle between J and the line of sight. The semi-major axis of the 
pulsar is defined as 

x^^^. (62) 



c 

where Op is the semi-major axis of the pulsar. By Eq.(59), we have, 

Op cos idi A • ^ • 

xi = — ;- = — xiismALjsmSicotz . (63) 

c at 

The semi-major axis of the orbit is a = ^cip, and since the L-S coupling induced d is a function 
of ft and uj, as shown in the appendix, we have 

d„sini a 

±2 = — = -X . (64) 

c a 

Therefore, the L-S coupling induced x is given by 



i: = xi -|- i;2 = — xJisin Alj sinOcot i -I — x . (65) 



By Eq.(65) we have 



X = xi I — + Q, cot H — — — - — — I -|- x- — — h X2- • (66) 
\n smAij J a'' a 

Notice that O, and O, can be obtained by Eq.(49). Considering Alj ^ 1 and by Eq.(60), Eq.(61), 
the observational advance of precession of periastron is given (Smarr &; Blandford 1976, Wex & 
Kopeikin 1999), 

a;°bs = u; + n- XLjcotisinn . (67) 

Therefore, we have 

cj°bs = ijj + n- Xl J cot i Sinn . (68) 

If UJ and $7 are caused only by the S-L coupling effect {H = Hg), then d> is given by Eq.(57). If we 
consider all terms of the Hamiltonian, as given by Eq.(2), then u> should include ci;*^^, the advance 
of precession of periastron predicted by general relativity, which caused by Hipjsf and H2pn- In 
such case to in Eq.(68) is replaced by w*^^ -|- to. Thus Eq.(68) can be written as 

^obs ^ ^GR ^ _ (gg) 

where 

uj^ = u! + tl- Xlj cot i sin n . (70) 

Notice that uj^ is a function of time due to Co and Xlj are function of time, as shown in Eq.(50), 
Eq.(57) and Eq.(99) respectively. Whereas, is a constant as shown in Eq.(8). 
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For a binary pulsar system with negligibly small eccentricity, the effect of the variation in the 
advance of periastron, co, is absorbed by the redefinition of the orbital frequency. As discussed by 
Kopeikin (1996), lo°^^ + Ae{u) IS given 

u;°^- + Ae{u)=uo + ^{t-to) , (71) 

where Ae{u) is the true anomaly, related to the eccentric anomaly, u, by the transcendental equation, 
loq is the orbital phase at the initial epoch Iq. 

At the time interval, 5t = {t — to), there is a corresponding 5io°^^ which causes a corresponding 
dPb at the right hand side of Eq.(71), therefore, Pj, is a function of time. Thus we have 

5cj°''' + Aeiu) = '^St . (72) 

Write l/Pb in Taylor series, we have 

Pb Pbito) piito) ^ ' 

Considering Aeiu) = 2'K5t/Pb{to) and by Eq.(72), Eq.(73) obtains, 

^^obs^_^^^_ (74) 
"ft 

Since Co^ is a function of time, whereas uj^^ =const, then we have iib°^^ = uj^ by Eq.(69). Assume 
F = tj'^, and write it in Taylor series as: F = Fq + F6t + ^FSt'^, obtains 6F = doj'^ ~ F5t = u^5t. 
Therefore, Sib°^^ of Eq.(74) becomes 6lo^ = io'^dt, from which Eq.(74) can be written as 

n = -^. (75) 

By Eq.(75), the derivatives of Pb can be obtained, 

p_2Pi Pid^.^^ P^d^uS 



d^P, _ Pi d^g;^ 
dt^ ^ 27r df^ ' 

6. Comparison of three different S-L coupling induced u"^^ and Pb 
6.1. Discrepancy between Wex &: Kopeikin and this paper 



(77) 



Sect 4 calculates the S-L coupling induced change of orbital elements of a binary system directly 
in the J-coordinate system (in which z axis is along J and x — y plane is the invariance plane). 
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And Sect 5 transform the effects in J-coordinate system to the observer's coordinate system, and 
obtains x, lo°^'^ and Pi,. In which a;"''* is equivalent to the definition of Wex &; Kopeikin (1999) as 
shown in Fig 2a. 

By Eq(57) the uj can be 1.5 PN (or w = in IPN), thus by Eq(70) the S-L couphng induced 
precession of periastron, io^ ^ becomes 

Cj^ tl — \l J cot ism^ . (78) 

Eq(78) is the result corresponding to 8 degrees of freedom, and uj^ is IPN. On the other hand, Wex 
& Kopeikin (1999), rewrite BO's orbital precession velocity, Eq(12) in the J-coordinate system, 
and then obtain i, ci;"'"' in observer's coordinate system. 

However the difference is that for Wex &: Kopeikin (1999), all the results in the J-coordinate 
system is calculated in to 10 degrees of freedom. In such case J in Eq(57) is replaced by S, thus the 
first term at the right hand side of Eq(57) is 0.5PN smaller than that of the second term. Therefore 
Co can be represented by the second term at the right hand side of Eq(57), which is IPN. By Eq(70), 
and considering tl — CIcosXlj ~ l.hPN {Xlj <S 1), we have 

Lo^ —Xi^j cot i sin U . (79) 

Although Xlj ~ S/L <^ 1, Xlj is significant which is IPN {Xlj ~ as shown in Eq(98) and 
Eq(99). Consequently lo'^ is also IPN, which leads to significant io^ (IPN), and therefore, significant 
derivative of by Eq(75)-Eq(77). In other words Wex & Kopeikin's u of Eq(59) (Wex & Kopeikin 
1999) actually corresponding significant variabilities, such as and Pf,, these results seem ignored. 

Therefore, this paper and Wex Kopeikin (1999), which corresponding to Eq(78) and Eq(79) 
respectively, both predict to significant ui^ and Pi,. However there is an obvious discrepancy between 
them, that is in Eq(79), uj^ when i 7r/2. Whereas, Eq(78) doesn't has such a relation. 
Therefore, the validity of Wex & Kopeikin (1999) and this paper can be tested by binary pulsar 
systems with orbital inclination, i — > 7r/2 (edge on). 

If a binary pulsar system with i it/2 still has significant P^ (IPN), then Eq(79) corresponding 
to Wex & Kopeikin (1999) is not supported, other wise Eq(78) corresponding to this paper is not 
supported. 

The discrepancy between Eq(78) and Eq(79) is due to the discrepancy on lj, which is caused 
by different degrees of freedom used in this paper and Wex & Kopeikin (1999). 

Eq(78) and Eq(79) have important property in common, that is uj'^ (and therefore, u>°^^) is 
obtained by the transformation from the J-coordinate system to the observer's coordinate system, 
in which all the three triads are at rest to SSB. 
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6.2. Discrepancy between Damour &; Schafer and Wex & Kopeikin 

Damour & Schafer (1988) express the orbital precession velocity as, 

where Kq unit vector along is line of sight, which defines the third vector of a reference triad (Iq, 
Jo, Kq), where Iq-Jq corresponds to plane of the sky. And the triad of the orbit is (i, j, k), in 
which k corresponds to L. i is the nodal vector determined by the intersection of the two planes 
(notice that it is different from the scaler, z, which represents the orbital inclination), as shown in 
Fig 2b. By Eq(80), and the relations between the reference triad, components of tig axe obtained 
(Damour k Schafer 1988). 

^ = ^[05-k-f25-Kocosi] , (81) 
at sin I 

| = n..i. (83) 

The S-L coupling induced Co given by Eq(81) is cj ~ Vts ~1-5PN. Therefore, Damour & Schafer 
(1988) predict insignificant to {u)^) which is much smaller than IPN, and therefore, the correspond- 
ing P{, is also insignificant. 

Thus it seems strange that Damour & Schafer (1988) and Wex &; Kopeikin (1999) start from 
same orbital precession velocity given by Eq(12), but predict different observational effect. 

This is because io^ of Wex &l Kopeikin is calculated in a coordinate system with axes {x' , y\ 
z') at rest to SSB as shown in Fig 2a. Whereas, io (w*^) of Damour k, Schafer is calculated in a 
coordinate system with triad (Kq, k, i), which is not at rest to SSB as shown in Fig 2b. Obviously i 
(point A), which is the intersection of the plane of the sky and the orbital plane of a binary system, 
is not static in the coordinate system (x', y', z') of Wex & Kopeikin (1999), so is k. In other words, 
triad (Kq, k, i) has non-zero acceleration relative to SSB, therefore, effects calculated based on 
such triad cannot be compared with observation directly. 

Obviously, if Damour & Schafer's Co is also calculated in the coordinate system as that of Wex 
h Kopeikin, then Eq(81) reduces to Eq(79). Thus, the discrepancy between Damour k Schafer 
(1988) and Wex k Kopeikin (1999) is that the former calculated in a coordinate systems which is 
not at rest to SSB; whereas, the latter is at rest to SSB. On the other hand, the discrepancy between 
Wex k Kopeikin (1999) and this paper is different degrees of freedom used in the calculation of 
equation of motion of a binary system. The relationship of the three kinds of S-L coupling effects 
is shown in Fig 4. 
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7. Confrontation with observation 

The precise timing measurement on two typical binary pulsars, PSR J2051-0827 and PSR J1713+0747, 
provides evidence on whether iu^ and Pj, is IPN or 1.5PN, but it still difficult to distinguish which 
IPN effect is valid, Wex & Kopeikin (1999) or this paper. 

The orbital motion causes a delay of T = ri • Kq/c = ri(t) sina;°^®(t) sini(t)/c in the pulse 
arrival time, where r is the pulsar position vector and Kq is the unit vector of the line of sight. 
The residual ST = ri ■ Kq/c — (ri • Ko/c)k of the time delay compared with the Keplerian value is 
of interest (Lai et al. (1995)). Averaging over one orbit r a, and in the case t <ti 1/\uj°^^\, the 
S-L coupling induced residual is, 

ST = ^ cos4tsina;°b« + ^E^tsmu"'°' + ^E^u^^H cos uj°'>' 
c at c c 

= xit sin 0;°''^ + -xt sin + u^'^'xt cos . (84) 
a 

By Eq.(94), we have a/ a ~ J/a^ ~ oj'^^^, therefore, the second and third term at the right hand 
side of Eq.(84) cannot be distinguished in the current treatment of pulsar timing. In other words, 
the effect of a can be absorbed by uj°^^. 



7.1. PSR J2051-0827 

As discussed above, the second term at the right hand side of Eq.(65) can be absorbed by uj°^^, 
therefore x ^ xi, and by Eq.(63), we have 

• X tani di 1 

\ u — — . (oO) 

X sin A Lj sin Qq dt sin A lj sin 

According to optical observations, the system is likely to be moderately inclined with an inclination 
angle i ~ 40° (Stappers et al. 2001). By the measured results of a; = 0.045 s, x = -23(3) x 10^^^ 
(Doroshenko et al. 2001), and by assuming sinA^jsinJlo = 2 x 10~^, Eq.(85) can be written in 
magnitude, 

n - ( ^ )( ^ )-V ^^^^ sinA^jsinOo , , , 

^^"4.3x10-13^^0.045^ ^tan40°^^ 2x10-3 ' ^ x lU (s j . i«bj 



In the following estimation of this section all values are absolute values. By Eq.(57), we can assume 
cj^ ~ (uj^)^ ~ 4 X 10^1^. Usually ui^ can vary in a large range, i.e., ij^ > {u!^)"^, depending 

on different combination of parameters, such as binary parameters, magnitude and orientation of 
5*1 and 5*2. In this paper we assume that io'^ ~ i^^)"^- Thus from Eq.(75) we have. 
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By Eq.(70), we can estimate d^u^/dt^ ~ Ri 8 x 10 s ^, similarly, we can estimate d^u^/dt^ ~ 
$7^ 16 X 10-36 

Therefore, by Eq.(76) and Eq.(77) we have A ~ 9 x IQ-^Og-i and d^Pb/dt^ ~ 2 x IQ-^^s-^. 
By Eq.(65) and Eq.(66), x/x ~ O ~ 2 x IQ-^s"^, which is consistent with observation as shown in 
Table 2. 

Therefore, once x is in agreement with the observation, the corresponding (jj^ can make the 
derivatives of be consistent with observation as shown in Table 2. Whereas, the effect derived 
from Damour & Schafer's equation cannot explain the significant derivatives of P^. 



7.2. PSR J1713+0747 

By the measured parameters, x = 32.3s, \x\ = 5(12) x 10"^^, i = 70° (Camilo et al. 1994), 
and by assuming sinA^jsinilo = 1 x 10"^, then in magnitude we have, 

^5 X 10-15 ''^32.3^ Han 70°^^ 1x10"^ ^ V y > \ J 

similarly we have. 

The comparison of observational and predicted variabilities are shown Table 3, which are well 
consistent. Notice that and P^*^** measured in these two typical binary pulsars cannot be 
interpreted by the gravitational radiation induced x and P^, since they are 3 or 4 order of magnitude 
smaller than those of the observational ones. 



7.3. PSRs J0737-3039A and B 

PSRs J0737-3039A and B is a double-pulsar system with P^ = 0.102251563(l)day, advance of 
periastron, cb = 16.90(l)deg/yr and orbital inclination angle i = 87.7I29 (^^rgay et al. 2003, Lyne 
et al. 2004). This binary pulsar system with z 7r/2 may tell us not only whether and P^ is 
significant or not, but also which significant 0;*^ and Pj, effects (Wex & Kopeikin or this paper) is 
valid. 

As given by Eq.(69) the measured lo is the sum of relativistic advance of periastron and the 
S-L coupling induced advance of periastron, 

cj°bs ^ ^GR ^^S ^ -^Q QQ (deg/yr) . (90) 
d)'^^ and io^ are both IPN. In order of magnitude one can estimate 



^ uj°^^ = 16.90 {deg/yr) 



(91) 
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Thus we have ~ {uj^)'^ 8.8 x lO^^^s-^. By the same treatment of PSR J2051-0827 and 
PSR J1713+0747 we have, 

The treatment ~ lo°^^ and lu^ ~ i^^)^ might over estimate the S-L coupUng induced Pb for one 

or even two order of magnitude. Nevertheless, the S-L couphng induced Pb is Hkely much larger 
than that caused by the gravitational radiation, P^^^ = —1.2 x 10~^^ ss~^ (Burgay et al. 2003). 
Since the observational Pb will be given soon(Burgay et al. 2003), whether Pb is significant or not 
can be tested on this binary pulsar system. 

The specific orbital inclination of this binary pulsar system can tell us more about the S-L 
couphng induced effects. By i = 87.7+29 of PSRs J0737-3039A and B, we have cot i ~ 0.04, therefore 
Wex and Kopeikin's cj^ given by Eq.(79) should be much smaller than that of this paper's lo^ given 
by Eq.(78). 

Thus we have two ways to test the validity the S-L coupling effect given by Wex and Kopeikin 
(1999) and this paper. 

The first one is that by u)^ ~ i'^^)^ oc (cot i)^, and Pb cx (w'^)^, we have {Pb)wK ~ (cot i)^Pb- 
Thus the magnitude of Pb corresponding to Wex and Kopeikin should not exceed 1.6 x 10~^^ ss~^, 
which is close to gravitational wave induced Pb- 

In other words if the measured magnitude of Pb of PSRs J0737-3039A and B is not much larger 
than that gravitational wave induced Pb, then the expression of this paper can be excluded. 

The second one is based on Eq.(75) to Eq.(77), from which we have, 

\Pb\ \Pb\ 1 .Qo^ 

\d^Pb/dt^\ ^ \P,\ ^ ^^^^ 

Since lo'^ corresponds to Wex and Kopeikin is nearly two order of magnitude smaller than the 
relativistic advance of periastron, lo'-'^; whereas, lj^ of this paper is same order of magnitude of 
io'-^^. Therefore, if the ratio of Eq.(93) is close to I/lo'-^^ then this paper is supported. 

The precise measurement of derivatives of Pb of this binary pulsar system may provide chance 
to test which one is favored. 



8. Discussion and conclusion 

BO and AK's orbital precession velocity has been treated as equivalent, since BO and AK give 
equivalent torque, L, as shown in Eq.(9) and Eq.(18) respectively. However Eq.(9) and Eq.(18) 
actually correspond to two different orbital precession velocities, Eq.(12) and Eq.(19) respectively, 
this is because the same torque can cause different effect when a dynamic system is calculated 
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under different degrees of freedom. BO's orbital precession velocity is actually obtained under 10 
degrees of freedom; whereas, AK's is under 8 degrees of freedom. Correspondingly the former one 
violates the triangle constraint and the latter one satisfies it. 

The discrepancy in physics results discrepancy in observation. Eq.(12) and Eq.(19) correspond 
to different combinations of 0, and to (fi and lu are defined in Fig. 2), and since observational effect 
depends on Q and to instead of L, therefore, the equivalent value in L may correspond to different 
observational effects. Concretely, BO and AK give same but different co (notice that Cl and co are 
components of the vectors given by Eq.(12) or Eq.(19)). And by Eq.(67) and Eq.(68), the observed 
advance of precession of periastron depends on both J7 and uj, thus BO and AK must correspond 
to different observational effects. 

In the calculation of Sect 4 and Sect 5, we can see the influence of the degree of freedom and 
physical constraint on the results of equation of motion, perturbation and observational effects. 

The S-L coupling induced precession of orbit can cause an additional time delay to the time of 
arrivals (TOAs), which can be absorbed by the orbital period. And since the additional time delay 
itself is a function of time, therefore orbital period change, P5, appears. Actually Pb corresponds 
to u)^ as shown in Eq.(75), which cannot be absorbed by cb'-'^ (d)^ can be absorbed by cu'-'^). 
Therefore, the higher order of derivatives of orbital period provide good chance to test different 
models. The observation of Pb, Pb and d^Pb/dt^ on PSR J 2051-0827 supports significant S-L 
coupling induced effects. 

This paper for the first time points out that Wex &; Kopeikin's expression in 1999 actually 
corresponds to significant cj^ and Pb, however it is not equivalent to the significant 0;*^ and Pb given 
by this paper. Precise measurement of P5, Pb and d^Pb/dt^ of specific binary pulsars with orbital 
inclination i vr/2, like PSRs 0737-3039A and B, may provide chance to test the validity of the 
results corresponding to Wex k. Kopeikin (1999) and that of this paper. 

I thank T. Huang for help in clarifying the theoretical part of this paper. I thank R.N. Manch- 
ester for his help in understanding pulsar timing measurement. I thank T. Lu for useful comments 
during this work. I thank W.T. Ni and CM. Xu for useful suggestions in the presentation of this 
paper. I thank E.K. Hu, A. Riidiger, K.S. Cheng, N.S. Zhong, and Z.G. Dai for continuous encour- 
agement and help. I also thank Y. Li, Z.X. Yu, CM. Zhang, L. Zhang, Z. Li, H. Zhang, S.Y. Liu, 
X.N. Lou, X.S. Wan for useful discussions. 

9. Appendix 

By 5 = a^o • n and T = a^o • t, ^ and ^ have been given by Eq.(49), Eq.(50), Eq.(55), 
and Eq.(57), following the standard procedure for computing perturbations of orbital elements 
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Roy ((1991)). Similarly, four other elements can be given: 

§ = ^i^(5esin/ + ^) , (94) 
at nv 1 - r 

da 4aiJ(l + e^) _ 
^ dt ^~ (l-e2)5/2a2^^^^^ ^^^^^ ' ^^^> 

de ^ VI -e2 ~ ^.^ ^ ^ ~ ^ ^ ^ 
at na 

<l>=(r^^W^''^^^"''^^^^' ^''^ 

c?Alj ^ I^rcos(a; + /) 1 
dt na^Vl - smXij ' 

at 2a-^[l — e^j-'/^ sm A^j 

+ {P,Qx - Pa,Qz){Syai + S2y<72) + (PxQy - PyQx){S,ai + ^2,(72)] , (99) 



^ = + 2^(1 - e2)^/^sin2 ¥) - ^ ' (100) 

dt 1 + vT^r^ dt dt^ ' ^ 2 ' na? ' ^ ' 

where ^ = ^ + 2f (sin^Y). 

< ^ >= ^ <^>+2<^>(l- e^)V2(3i,2 

dt 1 + ^1 _ e2 dt dt ^ ^ ^ 2 ^ 

4(71 J 



a3(l - e2 



■(P^g^/ - -PyQa:) • (101) 



REFERENCES 

Apostolatos,T.A., Cutler, C, Sussman,J.J. Thorne,K.S., 1994, Phys. Rev. D, 49, 6274-6297 . 
Barker,B.M., O'Connell, R.F., 1975, Phys. Rev. D, 12, 329-335. 
Burgay, M. et al. 2003, Nature, 426, 531. Astrophys. J. 437, L39-L42 . 
Camilo, F., Foster, R.S., Wolszczan, A., 1994, Astrophys. J. 437, L39-L42 . 
Damour, T., Schafer, G., 1988 IL Nuovo Cimento, lOlB, 127 . 

Doroshenko,0., Lohmer, O., Kramer, M., Jessner, A., Wielebinski,R., Lyne,A.G., Lange, Ch., 2001, 
Astro-Astrophys, 379, 579-588 . 



-23- 



Hamilton, A.J.S. Sarazin, C.L., 1982, MNRAS 198, 59-70 . 

Kidder, L.E., 1995, Phys. Rev. D, 52, 821-847. 

Kopeikin, S.M., 1996, Astrophys. J. 467, L93-L95. 

Lai, D., Bildsten, L., Kaspi,V., 1995, Astrophys. J. 452, 819-824 . 

Liu, L., 1993, Method in Celestial Mechanics, Nanjing University Press 

Lyne, A. G. et al. 2004, Science, 303, 1153. 

Roy, A.E., 1991, Orbital motion, Adam Hilger 

Smarr, L.L.,Blandford, R.D., 1976, Astrophys. J. 207, 574-588. 

Stappcrs, B.W. , Van Kerkwijk,M.H., Bell, J.F., Kulkarni, S.R., 2001, Astrophys. J. 548, L183- 
L186 . 

Wex, N., Kopeikin, S.M., 1999, Astrophys. J. 514, 388. 

Yi, Z.H., 1993, Essential Celestial Mechanics, Nanjing University Press 



This preprint was prepared with the AAS IM^^X macros v5.2. 



-24- 



Table 1: Comparison of S-L coupling induced variabilities given by different authors 

PS (1988) WK (1999) This paper evidence 

n in J-co IPN IPN 

u in J-co IPN 1.5PN 

oj^ ~ (of Eq.(12)) — Aijcotisinn Cl — Xlj cot ismU 

^obs a;G^+1.5PN cL-^^+lPN l:;^^+1PN 

Pf) pGR l^obsj ^ ji^'-'^l l^obsj ^ j^G-Rj |^obs| ^ 

wheni^7r/2 pobs ^ pGR jpofe^j ^ |pGR| 

J-co represents J-coordinate system, P^^ represents orbital period change due to gravitational radiation 
predicted by General Relativity, tls is given by Eq.(12) which is 1.5PN, whereas, tl and \lj are IPN, as 
given by Eq.(50) and Eq.(98) respectively. 



Table 2: Measured parameters compared with the geodetic precession induced ones in 
PSR J2051-082 7 

observation WK & this paper 



f» -23(3) X 10-1^ X = x°''" 

{x/xf^'' < -3.0 X lO^^s"^ \x/x\ PS 2 X lO-^s"^ 

P°bs = -15.5(8) X 10^12 |p^| ^ jii^i _ 5 ^ ;^Q-ii 

P^^' = 2.1(3) X 10-20s-i = l^^^l ~ 9 X 10-20s-i 

^ = 3.6(6) X 10-^^s-^ 1^1 = fg^l ^ 2 X 10-^^s-^ 



Table 3: Measured parameters compared with the geodetic precession induced ones in 
PSR J1713+0747 ^^^^^^^^^^^^^^^^^^^^^ 

observation WK & this paper 

|x|°'^^ = 5(12) X 10-1^ X = 

p°bs ^ ^29) X 10-1^ inl = \^^\ ~ 1 X 10-1° 
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Fig. 1. — In IPN, the scenario of motion of a binary pulsar system can be described as the rotation 
of the plane determined by L, S around the fixed axis, J. ry^ and rjs are precession phases of L 
and S in the J-coordinate system respectively. 





Fig. 2. — Fig 2a: binary geometry and definitions of angles of Wex &: Kopeikin (1999) and this 
paper. The invariable plane (x-y), as represented by the dotted ellipse, is perpendicular to the 
total angular momentum, J. The inclination of the orbital plane with respect to the invariable 
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Fig. 3. — Angles and orientation conventions relating vectors S, L and J to the coordinate system, 
x-y is the invariance plane, risi, ris2, Vs and rjL are precession phases of Si, S2, S and L, respectively. 
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